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ABSTRACT 

Local and average Nusselt numbers have been experimentally 
determined around the circumference of a uniformly heated 
right circular cylinder placed near a plane surface in a 
cross flow of air. Experiments were carried out for Reynolds 
numbers ranging from 16,000 to 153 . 000 , three plate lengths 
and varying cylinder to plate spacings. The variations of the 
average Nusselt number as a function of gap spacing was found 
to exhibit an interesting and unexpected trend. The average 
Nusselt number exhibited a minimum value when the cylinder 
was resting on the plate and attained a maximum value at a gap 
spacing of approximately one to two cylinder radii. Flow 
visualization studies together with measurements of the vortex 
shedding frequency in the cylinder wake indicate that when the 
plate is within approximately one cylinder radius of the 
cylinder, it inhibits the formation of the regular vortex 
pattern that forms in the wake of a free cylinder . This has 
a strong influence on the base Nusselt number which, in turn, 
exerts a controlling influence on the magnitude of the average 
Nusselt number. 
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I. INTRODUCTION 



An experimental investigation was carried out to examine 
the influence of wall proximity effects on the heat transfer 
characteristics of a uniformly heated cylinder placed in a cross- 
flow of air. Recently, McComas /~ 1_J7 initiated an investigation 
of the heat transfer characteristics of a uniformly heated 
cylinder attached to a plate and Gnerlich 2_/ continued the 
same study varying the cylinder to plate spacings. The present 
work is an extensive follow-on study of the work of McComas 
and Gnerlich differing in that a range of Reynolds numbers has 
been investigated as well as three different plate lengths. 

The present experimental work was conducted in the Mechanical 
Engineering Department wind tunnel (described in detail in 
Section III) in Building 500 at the Naval Postgraduate School. 

Local Nusselt numbers were obtained around the circumference 
of the cylinder at Reynolds numbers of 16,000, 90,000 and 
153,000 without a plate as well as with a plate. In the 
case with a plate, the cylinder to plate spacings were varied 
from d/r = 0 to d/r =5-33 with L/D = 4 where d = gap width, 
r = cylinder radius, D = cylinder diameter and L = length of 
plate from the centerline of the cylinder to the leading edge 
of the plate (Figure 1) . Experiments were also conducted at 
a Reynolds number of 153,000 with L/D = 2 and 8 in order to 
see the effect of the length of the plate's leading edge on 
the heat transfer results. 
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Figure 1. Cylinder-plate configuration. 



Before proceeding further, a "brief qualitative description 
of flow past a cylinder will help to explain and interpret the 
results obtained from this experimental study. Since the 
present experiments examined Reynolds numbers ranging from 
16,000 to 153,000, only the subcritical and critical flow 
regimes will be described here. See Achenbach /~3_7 for a 
complete discussion of the various flow regimes. 

In subcritical flow, a laminar hydrodynamic boundary 
layer grows from a minimum thickness at the forward stagnation 
point on the cylinder to a maximum thickness at an angular 
location of approximately 80-85 degrees. At this location, 
the kinetic energy of the moving fluid is no longer sufficient 
to overcome the adverse pressure gradient present on the 
cylinder’s surface and the laminar boundary layer separates. 
Upon separation a turbulent wake then develops on the rear 
side of the cylinder. 

In critical flow, the laminar boundary layer again grows 
from the forward stagnation point to the laminar separation 
point which occurs at approximately 80-85 degrees. In contrast 
to the subcritical flow case, a transition occurs from laminar 
to turbulent flow. A turbulent boundary layer then reattaches 
itself to the cylinder. This reattached turbulent layer 
eventually separates at a position downstream of the laminar 
separation point. The region between the laminar separation 
point and the point of reattachment of the turbulent boundary 
layer is known as a "separation bubble." 
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The results of the experiments with the cylinder attached 
to the flat plate indicate a critical flo w pattern exists 
on the top half of the cylinder for all Reynolds numbers 
investigated. This flow pattern changed from critical to 
subcritical over the top as the cylinder to plate spacing 
was increased to approximately d/r = 0.1 to 0.2. The flow 
over the bottom half of the cylinder revealed critical flow 
immediately after a space existed. Flow over the bottom 
turned from critical to subcritical in the interval of 
0.8 - d/r - 1.2. 

Plots of the average Nusselt number versus gap spacing 
have shown that the minimum average Nusselt number occurs when 
the plate is attached to the cylinder and the maximum average 
Nusselt number occurs at d/r = 2.0. 

Qualitative studies /~ 4 _ 7 , using a flow visualization 
technique in the water tunnel together with measurements of 
the vortex shedding frequency in the wind tunnel, have contri- 
buted a plausible explanation about the trends of average 
Nusselt number curves. 

It is theorized that trends in the average Nusselt 
number curves are a direct consequence of the strength and 
frequency of the vortex pattern which forms in the wake of 
the cylinder. 
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II. BACKGROUND 



The present investigation of the heat transfer character- 
istics of a uniformly heated cylinder placed in close proximity 
to a plane surface in an air stream was conducted concurrently 
with an investigation by Goktun /~4_7 of the flow character- 
istics of a cylinder placed near a plane surface. Information 
obtained from the survey of the heat transfer experiments 
proved useful in explaining the nature of the flow field 
observed by Goktun. In the same manner, Goktun' s observations 
made significant contributions to a better understanding of 
the heat transfer characteristics of the cylinder. 

Although extensive investigations have been directed at 
obtaining a better understanding of the heat transfer and flow 
characteristics of "free" cylinders (for extensive references 
on this subject see Cooper, Field and Meyer /T 5_7) < data 
relating to the heat transfer and flow characteristics of 
the cylinder-plate geometry under investigation were not 
available in the literature. Such information would be of 
obvious value in the design of heated piping systems placed 
near walls and in the design of heat exchangers. 

As previously mentioned, the present investigation is 
a follow-on study of an investigation initiated by McComas 
7~ 1_7 and continued by Gnerlich 2_7 • Although McComas 
experimentally determined local and average heat transfer 
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coefficients around a uniformly heated cylinder attached to 
a wall and Gnerlich did likewise for various cylinder to 
plate spacings, "both of these investigators found no references 
to theoretical work on this subject. 

Convective heat transfer is significantly affected by 
the fluid flow characteristics occurring adjacent to a solid 
boundary. Since forced-convection heat transfer is a result 
of the movement of fluids and the mixing of fluid particles, 
the mechanism of fluid flow must be understood in order to 
understand the mechanism of heat transfer and to explain the 
phenomena occurring during the heat transfer process. Because 
of this, this project was conducted concurrently with Goktun's 
/~4_7 study of the flow field near the cylinder-plate geometry. 
In explaining the heat transfer results obtained in this 
work, frequent reference will be made to Goktun's concurrent 
flow investigation. 
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III. EXPERIMENTAL APPARATUS 



1. Wind Tunnel 

An open circuit, subsonic wind tunnel was used for the 
experiments. The tunnel was specially designed by Kenney 
Engineering Corp., of Monrovia, California for research 
quality work. The entire exterior dimensions are given in 
Figure 2. The tunnel has a test section which is 20 inches 
wide, 28 inches high at the entrance and expands to 31 inches 
at the exit to provide for boundary layer growth. It is 
powered by a 75 HP, 1?50 RPM variable speed controlled D.C. 
motor which provides a belt drive to an adjustable pitch 
vaneaxial fan. The fan is 45 inches in diameter and is 
capable of supplying air at 10,000 to 70,000 cfm, producing 
wind velocities ranging from 10 mph (14.7 ft/sec) to 175 mph 
(256.7 ft/sec). The tunnel is equipped with six graded mesh 
screens to control turbulence. A 24 inch water manometer 
connected in to the stilling chamber and test section may be 
used to measure velocity. Utilizing a hot wire anemometer, 
a turbulence intensity level of 0.2 % was found in the wind 
tunnel at a free stream velocity of 98 . 1 ft/sec. A mercury 
thermometer was installed in the stilling chamber to measure 
the free stream temperature. 
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Figure 2. Schematic diagram of the wind tunnel. 




2. Liquid Crystals 

Microencapsulated cholesteric liquid crystals were used to 
obtain both qualitative and quantitative heat transfer and 
fluid flow information on the uniformly heated cylinder. The 
liquid crystals provide a visual display of temperature over 
discrete, reproducible temperature ranges (called the "event" 
temperature range) by means of a change of color (see Cooper, 
Field and Meyer 5 _ 7 ) • The visible color spectrum, which 
falls within a certain temperature band, has an accuracy of 
0.2°F. All experiments in this investigation were conducted 
using a liquid crystal designated as NCR S-^3 having an event 
temperature range of 

Red : 109. 8°F ( 43 . 2°C ) 

Green: 111.0°F (43„9°C) 

Blue : 112. 3°F (44.6°C) 

3 . Temsheet 

An electrically resistive carbon impregnated paper known 
as "Temsheet" was used as the heating element in these studies. 
Temsheet is a thin, highly flexible, fibrous material contain- 
ing no wires or ribbons and therefore is free of localized hot 
spots, burnout, and breakage problems. The nominal thickness 
is 0.1 cm and the electrical resistivity is about 2.5 ohm-cm. 
Using the Temsheet it is possible to heat large areas within 
less than 2 °/o deviation from point to point. 
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4. Cylinder 

The test cylinder utilized in obtaining the heat transfer 
data was the one used by Gnerlich 2_7. It was made of three 
cylinders attached end to end (see Figure 3) having a total 
length of 20 inches, corresponding to the test section width, 
and a diameter of 3 inches. The two end sections were hollow 
wood cylinders connected by a wooden brace, curved to match 
the cylinder contour. The central cylinder, henceforth 
designated as the test section, was formed from a sheet of 
Temsheet. It was tightly packed with glass wool which served 
the dual purpose of providing insulation against possible 
internal heat losses due to conduction and free convection 
and also provided the paper cylinder with added resistance 
to deformation. The circumference of the Temsheet cylinder 
was marked at five degree angular intervals in two locations 
2 inches apart for accurate location of the liquid crystal 
isotherms during the experiment. Aluminum foil conducting 
tape served as electrodes. The aluminum tape was circumfer- 
entially attached on the inner surface of the Temsheet. Silver 
conducting paint was applied to the interface between the tape 
and Temsheet to insure good electrical contact was maintained 
between the two surfaces. Two lead wires were connected to 
the aluminum foil through one end of the cylinder. The lead 
wires were attached to a regulated DC power supply (model: 
LAMBDA LK 345A FM) . The resistance of the Temsheet was 
measured by a digital multimeter (model: FLUKE 8100A) . 
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Figure 3 Photograph of the Temsheet cylinder. 



5. Plates 

The plates used in the experiments were constructed of 
plexiglass and were 0.5 inch thick, 20 inches wide and had 
total lengths of 40, 46 and 58 inches, respectively. The 
plates were installed horizontally in the wind tunnel in such 
a way that the leading edge of the plates extended 2, 4 and 
8 cylinder diameters (i.e. 6, 12 and 24 inches, respectively) 
ahead of the cylinder centerline. The plates were beveled 
on the under side of the leading edge with an angle of 12 
degrees. The plates were capable of being located at any 
position relative to the cylinder (Figure 4), from attached to 
a maximum distance of 8 inches away (i.e. d/r = 5*33)- The 
plates were moved vertically by means of a threaded rod 
extending through the tunnel floor, and a rear brace which 
locked in position by means of a split nut. 
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Figure 4 Photograph of the plate and cylinder assembly. 






IV . EXPERIMENTAL PROCEDURE 



The present investigation was conducted using the same 
test facilities as employed by Gnerlich 2_7 • Gnerlich, 

however, performed experiments only for a Reynolds number of 
153.000 (based on the cylinder diameter) and a plate leading 
edge of 12 inches. In the present work, experiments were 
conducted at Reynolds numbers of 16,000, 90,000 as well as 
153.000 and the plate leading edge was varied from 6 to 24 
inches . 

To initiate the experiment, a sufficient voltage (from 
12 to 15 volts) was applied to the cylinder to insure that the 
liquid crystals reached their event temperature and changed 
color over the entire surface of the cylinder. This was done 
to check for hot and cold spots which would indicate improper 
packing of the glass wool inside the cylinder. The wind 
tunnel was then started and set to a predetermined velocity 
at which it was maintained throughout the run. The applied 
voltage was adjusted until the first faint red color trace 
was observed. The power supply was then turned off and the 
multimeter was switched from DC voltmeter to ohmmeter. The 
same procedure was also followed at the end of the experiment. 
In this manner, an accurate value of the paper resistance, 
which was taken to be the arithmetic mean of initial and 
final resistance values, was determined. The multimeter was 
then switched back to the DC voltmeter and power was re-applied 
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to the test section. A steady state condition was usually 
attained in from five to ten minutes depending upon the 
wind velocity. Generally, at least ^-0 minutes were spent 
to take the first data point in order to be extremely certain 
about steady state conditions. 

The first color band observed on the cylinder was at the 
location of the value of the lowest local Nusselt number. 

The voltage was then increased such that the liquid crystal 
color bands would appear at various angular positions around 
the circumference of the cylinder. The temperature obtained 
in this manner could be compared to the air stream temperature 
which was obtained from a thermometer mounted in the stilling 
chamber. Using this temperature difference and power dissi- 
pation, the local Nusselt number based on the cylinder diameter 
could be obtained around the cylinder (see Appendix A). 

Due to the method of construction, there was a seam along 
the length of the test cylinder. Measurements within 30 
degrees of a sector of the seam were considered unreliable. 

The majority of Nusselt number data were collected with the 
seam located at 180 degrees. The Nusselt number data in the 
vicinity of .180 degrees were taken by turning the cylinder 
such that the seam was facing forward at 0 degrees. Several 
other seam positions were also checked but it was ultimately 
decided that the 0 and 180 degrees locations were the best 
for the seam position. In all experiments, the location of 
the liquid crystal color transition points were easily 
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distinguishable around the circumference of cylinder except 
in the region of 0 to 45 degrees and 160 to 200 degrees. 
Especially in the forward stagnation region the temperature 
gradient was very shallow and the whole region appeared 
isothermal. In the vicinity of the rear stagnation point, 
the temperature gradient was indistinct due to the turbulent 
behavior of the flow in the cylinder wake . 
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V . RESULTS 



The results shown in this section were obtained from 
the heat transfer experiments discussed in Section VI . The 
results are classified in two groups: 

A. LOCAL NUSSELT NUMBERS, Nu 

a. Typical data and a sample calculation for the heat 
transfer coefficient and the local Nusselt number are shown 
in Appendix A . 

b. Figures 5 through 12 show the variation of the local 
Nusselt number with gap size at a Reynolds number of 153 » 000 
and a plate leading length of 12 inches (L/D = 4 ) . These 
plots also show the variation of the pressure coefficient 

as a function of gap size. All pressure data we re taken 
from Goktun's /~4_7 results and are informative in that they 
show the relationship between surface pressure and heat 
transfer characteristics. 

c. Additional plots of the Nusselt number versus gap 
size for Reynolds numbers of 16,000, 90,000 (with L/D = 4 ) 
and 153 » 000 (with L/D = 2 and 8 ) can be found in Appendix D, 
Figures 24 through 55- Local Nusselt number distributions on 
the surface of a free cylinder are also contained in Appendix D, 
Figures 56 through 58. 

B. AVERAGE NUSSELT NUMBERS, Nu 

a. A graph of TTu versus d/r (the cylinder to plate 
spacings) for Reynolds numbers of 16,000, 90)000 and 153 >000 
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and L/D = 4.0 is shown in Figure 13 . A comparison is made 
with Gnerlich's £~ 2_7 results on the same graph. 

b. A graph of Nu versus d/r for Reynolds number of 
153,000 and L/D =2, 4 and 8 is shown in Figure 14. 

c. A graph of normalized drag coefficient (C,/C, ), 

Strouhal number (S/So, ) and average Nusselt number (Nu/Nuo, ) 
versus d/r is shown in Figure 15. The "co" subscript refers 
to values obtained with no plate in the tunnel. 

d. A graph of normalized Strouhal number and local 
base Nusselt number (Nu^^/Nu^ -b ase ) versus d/r is shown 
in Figure l6. 
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Figure 5- Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 
Re = 153,000, L/D = 4, d/r =0.0. 
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Figure 6. Local Nusselt number and pressure coefficient on 

the surface of a cylinder placed near a plane surface 
for 

Re = 15 000, L/D = 4, d/r = 0.25. 
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Figure 7 . Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 
Re = 153,000, L/D = 4, d/r =0.5- 
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Figure 8. Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 
Re = 153,000, L/D = 4, d/r =1.0. 
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Figure 9* Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 
Re = 153,000, L/D = 4, d/r = 2.0. 
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Figure 10. Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 
Re = 153,000, L/D = 4, d/r =3-0. 
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Figure 11 . Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface 
for Re = 153 , 000 , L/D = 4, d/r = 4.0. 
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Figure 12. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153,000, L/D = 4, d/r = 5-33. 
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length and gap size for a Reynolds number of 153 » 000 . 
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0.0 1.0 2.0 3.0 4.0 5.0 7.0 8.0 9-0 10.0 d/ / r 

Figure 15. Variation of normalized drag coefficient ) > Strouhal 

number ( S / ) and average Nusselt number (Nu / Nug, ) as 

a function of gap size for Re = 153 >000 and L/D = 4. 
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Figure 16. Variation of normalized Strouhal number ( S / Sqq ) and base 
Nusselt number (Nu base / Nu roi)age ) as function of gap size 
for Re = 153.000 and L/D = 4. 



YI . DISCUSSION 



Because of the relationship between the flow characteristics 
(i.e. pressure distribution around the cylinder, pressure drag 
and vortex dynamics) and the heat transfer characteristics of 
the cylinder, the parallel flow study of Goktun /" 4_7 will 
be referred to in some parts of the following discussion. 

Goktun obtained his results using an unheated aluminum cylinder 
placed in proximity to the plane surface. 

The general nature and trends of the results obtained 
from the present investigation can be discussed in two 
groups : 

A . TRENDS IN LOCAL NUSSELT NUMBERS , Nu 
At d/r = 0 

The initial heat transfer measurements were made with the 
uniformly heated cylinder attached to the plate at a Reynolds 
number of 153 )000 with L/D = 4. The data obtained in this 
arrangement agreed within the estimated experimental uncertainty 
of ±5# (See Appendix B) with the results of McComas /~1_/ 
and Gnerlich 2_J7 • A local maximum and minumum value of 
the Nusselt number were once again observed near the point 
of attachment of the plate and cylinder and are indicative of 
a large trapped vortex as proposed by McComas 1 _/ . In the 

flow visualization experiment, at d/r = 0, Goktun /”4_7 also 
observed that a series of small, circular, three dimensional 
vortices were also present in the upstream cavity as well as 
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a large turbulent vortex. The large turbulent vortex occurs 
within the separated region where the flow separates from the 
flat plate and reattaches to the cylinder. The point of flow 
reattachment occurs at approximately 140 to 1?0 degrees from 
the point of attachment of the plate and cylinder. The 
existence of this large turbulent vortex produces the observed 
local maximum and minimum values in the Nusselt number curve. 

The existence of a separation bubble on the top of the 
cylinder was the obvious mark of critical flow. A laminar 
boundary layer started to grow on the forward half of the 
cylinder and separated between 80-85 degrees and a subsequent 
transition to turbulent flow occurred resulting in the 
reattachment of a turbulent boundary layer at about 100°-105°» 
The reattached turbulent layer ultimately separated at 
approximately 130--135 degrees. 

The decrease in the local Nusselt number (see Figure 5) 
at the aft end of the cylinder was unexpected. The flow 
visualization studies clearly showed, however, that when the 
cylinder is resting on the plate a large wake with very little 
fluid motion develops on the downstream side of the cylinder. 
The absence of fluid motion therefore explains the decrease 
in Nusselt number on the back of the cylinder. 

At d/r = 0.23 to d/r = 5.33 

Further experimental results at a Reynolds number of 
153,000 were obtained with the plate as far removed from 
the cylinder as possible, d/r = 5-33* 
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Examination of the polar plots in Figure 6 to 12 for 
significant trends or unusual nature show two items of 
interest: 

1 . A shift in the location of the forward stagnation 
point, rear stagnation point and the laminar separation points 
( see Figure 17) . 

Starting from d/r = 0.25 and moving the plate to 
d/r = 2.0, it appears as though an imaginary axis initially 
passing through 340 and 160 degrees is pivoted about the 
origin. As this axis turns clockwise through 20 degrees of 
arc to 360 degrees, the angular position of the stagnation 
and separation points changed correspondingly up to d/r = 2, 
and then reached a constant behaviour (see Figure 6 to 9) . 

This observed shifting of stagnation points was 
definitely not affected by the variation of the length of 
the plate's leading edge and/or the Reynolds number as shown 
in Appendix D, in Figures 40 to 55 * 

2. The transition from critical to subcritical flow. 

The transition from critical to subcritical flow was 

investigated in order to find the d/r ratio associated with 
the transition on the top and bottom of the cylinder. At 
a Reynold's number of 153 ,000 and L/D = 4 , the flow was 
critical over the top of the cylinder with the plate attached 
and was subcritical over the top at d/r = 0.25> By alternately 
varying the plate to cylinder spacing, carefully adjusting the 
voltage and noting the liquid crystal display, the transition 
point from critical to subcritical flow was precisely found 
at d/r = 0.1. 
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Figure 17 . Shift of the Laminar separation points and the stagnation point 
as a function of d/r. 



A transition from critical to subcritical flow over 
the bottom of the cylinder was found at approximately 
d/r =1.2. In this case, although the pattern of the 
separation bubble displayed by the liquid crystal band was 
less pronounced than over the top, it was clearly observed 
as far as d/r = 1.2. 

As also noted by Gnerlich /~ 2 _7, flow over both the 
top and bottom of the cylinder are definitely critical for 
the stated d/r ranges but there is a significant difference 
in the degree of criticality on the top and bottom of the 
cylinder as far as the magnitudes of the local Nusselt numbers 
are concerned. This is seen in Figures 5 through 8. 

Further study was also conducted in order to see the 
influence of the length of the plate's leading edge and 
Reynolds numbers. The results are shown in Table 1. Table 1 
shows the transition points from critical to subcritical flow 
that occurred over the top and bottom of the cylinder as a 
function of Reynolds number and L/D ratios. It is worthwhile 
to note all points fell within the range of 0.1 - d/r - 0.2 
for the top. The transition points for the bottom were 
found in the interval of 0.8 - d/r - 1.2. Also to be noted 
in Table 1 is that the transition from critical to subcritical 
flow at various d/r ratios is apparently independent of 
L/D ratios but is mildly dependent on Reynolds number. As 
the Reynolds number increased, the d/r ratios for transition 
on the top of cylinder decrease slightly; on the bottom of 
the cylinder the behaviour was reversed. 
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0.15 
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0.1 
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Table 1 . Transition from critical to subcritical flow as a 
function of Reynolds number, d/r and L/D ratios. 

B. TRENDS IN AVERAGE NUSSELT NUMBER , Nu 

In order to gain further information about the cylinder 
heat transfer characteristics, the average Nusselt number, 

Nu was computed from the formula 

2 TT 

Nu = 2rf f Nu’de 
o 

The integration was carried out on an IBM-360 computer using 
Simpson's Rule of Quadrature with a step size of 5 degrees. 

Figure 13 shows the variation of the average Nusselt 
number as a function of gap spacing for three Reynolds 
number and L/D = 4 . 

A comparison of the results obtained at Re = 153 .000 and 
L/D = 4 , with those obtained by Gnerlich indicates that the 
trends are essentially the same. The decrease in Nu at 
d/r = 0.5 in Gnerlich' s results was not observed in the 
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present experiments but, considering the magnitude of the 
decrease and the 5^ experimental uncertainty involved in 
both experiments, it is not felt to be significant. The 
maximum Nu was observed at d/r = 2.0 and it is readily 
apparent that the presence of the plate radically altered 
the magnitude of Nu in the range 0 - d/r - 2. Beyond a 
d/r of 2.0, Nu approached a nearly constant value. 

It is important to emphasize that the trends of the Nu 
versus d/r curves are nearly the same for all Reynolds 
numbers . 

In order to see the effect of the length of the plate's 
leading edge, results were obtained for a Reynolds number of 
153 » 000 lor L/D ratios of 2, 4 and 8. The results, which are 
shown in Figure 14 , clearly indicated that there is no 
significant effect on the trends in the Nu curve due to 
changing the L/D ratios. Goktun /~ 4_7 also noted that 
changing the L/D ratio had no significant influence on the 
drag and lift characteristics of the cylinder. 

For the sake of simplicity, the flow characteristics 
(pressure drag and vortex shedding frequency) and the heat 
transfer characteristics (average and base local Nusselt 
number) are combined and plotted in normalized form in 
Figures 15 and 16 . 

Goktun' s qualitative flow visualization study /~ 4_7 in 
the water tunnel together with the vortex shedding frequency 
measurements in the wind tunnel indicated that vortex dynamics 
control the drag on the cylinder. Figure 15 clearly indicates 
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the interrelation between the drag and heat transfer. It is 
seen that as the drag coefficient increases, the average heat 
transfer coefficient also increases. It is therefore concluded 
that the vortex dynamics are responsible for the trends 
observed in both the drag and average Nusselt number curves. 

In the range 0<d/r<(l, Goktun found that the natural or 
regular vortex pattern formed in the wake of a free cylinder 
is inhibited by the presence of the plate. The plate signifi- 
cantly reduces the vortex shedding frequency. This in turn has 
an effect on the base pressure (see Figures 5 through 8) and 
base Nusselt number (Figure 16 ) . The more stagnant the wake, 
the higher the base pressure and the lower the base Nusselt 
number. Therefore, the stagnant wake results in poor heat 
transfer. This trend, in turn, influences the drag and 
average Nusselt number (Figure 15)* 

Goktun /“4_7 observed that strong and violent vortices 
were formed behind the cylinder when d/r was approximately 
1.0 to 2 . 0 . These vortices transfer high-momentum cool fluid 
from outside the boundary layer in toward the rear surface 
of the cylinder. The cool fluid in the vortex which impacts 
on the hot cylinder surface has a large capacity to remove 
energy. The "scrubbing" motion of these vortices building 
up behind the cylinder and breaking away to be carried down- 
stream in the wake, results in high rates of heat transfer 
and this results in a large base Nusselt number. This trend, 
in turn, influences both the drag coefficient and the average 
Nusselt number. 
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As the plate is moved further from the cylinder, the 
drag and Nusselt number approach their "free" cylinder 



values . 



VII. SUMMARY 



This experimental investigation was carried out to gain 
an understanding of the heat transfer behavior of a uniformly 
heated right circular cylinder placed near a plane surface. 

A concurrent study by Goktun /~A_7 was also carried out to 
study the lift and drag characteristics of a cylinder in the 
same configuration. Local Nusselt numbers were measured around 
the cylinder circumference for various plate to cylinder 
spacings, plate lengths and Reynolds numbers. The average 
Nusselt numbers were obtained from the local Nusselt number 
distributions by numerical integration. 

An unusual trend which was observed in the plots of both 
the average Nusselt number and overall drag versus cylinder to 
plate separation curves lead to a flow visualization experiment 
and hotwire anemometer measurements of the vortex shedding 
frequency in the cylinder wake. Based on these observations 
and measurements, it was concluded that the structure of 
vortex pattern that develops in the wake of the cylinder 
exerts a controlling influence on the shape of the average 
Nusselt number and drag curves. In the range of a plate to 
cylinder spacing of one cylinder radius, it was found that 
the presence of the plate inhibits the formation of the 
regular vortex pattern that forms in the wake of a free 
cylinder. This significantly affected both the cylinder 
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base pressure and local Nusselt number distribution which, 
in turn, influenced the drag and average Nusselt number 
curves . 
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APPENDIX A 



HEAT TRANSFER DATA REDUCTION 

The procedure followed in reducing the heat transfer 
data is shown below. In this manner the local Nusselt 
numbers at various locations were obtained for the case 
of 'a uniformly heated (constant heat flux) cylinder. 



HEAT FLUX 



The heat flux produced by Joulean heating is given by; 
^2 



% = 



A-R 



where : 



= Voltage applied to the test section, /~volts_7 
= Total resistance of Temsheet, /~ohms_7 

r\ 

= Total heated area of Temsheet ( =0.315) /""ft J 

.2 



w 



= Rate of total heat transfer, 



V‘ 



n- ft 4 



Conversion: 1 



V‘ 



a- ft' 



= 3-^13 



BTU 



hr-f t 



Hence , 



• _ 3.413 E' 

q w 0.315 R 



= 10.842 



R 



BTU 



hx'-ft' 
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HEAT TRANSFER COEFFICIENT 



The total heat transfer coefficient is determined by 
dividing the heat flux by the measured local temperature 
difference (T - T^ ) : 



h t = 



T - T 






CO 



AT 



h t = 



h + h 
c r 



h = 
c 



x w 



AT 



- h. 



h r = ' (T + T co ) (T 2 + T 2 ) 



where : 
h , 



T 



T« 




a 



total heat transfer coefficient, 

hr-ft 2 - °F 

cylinder surface temperature at angular 
location e , as determined by liquid crystals, 

air stream temperature, °F 



convective heat transfer coefficient, 
radiation heat transfer coefficient, 

Stef an-Boltzman constant = 0.1714 x 10 



BTU 

hr - ft 2 
BTU 

hr - ft 2 

-8 BTU 

hr - ft 




radiation exchange factor between cylinder 
and surroundings = 0.9 



NUSSELT NUMBER 



Nu = 



h • D 
c 



k 



f 



where : D 



cylinder diameter = ft 



k^ = thermal conductivity of air evaluated at 



NUMERICAL EXAMPLE: 

A sample calculation of the local Nusselt number is 
shown below using the data (seepage 57 ) obtained at a 

Reynolds number of 153 , 000 , d/r = 4.0, L/D = 8 and angular 
location of 6 = 63 degrees. 



V = 31 volts 
R = 13.3 

T = 111°F = 5?1 °R 
Tccp 70°F = 530°R 
iff = T - = 4l°F 



film temperature T^ 



BTU 



2 



’ hr - ft - °F 



3T-2 = °-9 

„ ,„-8 





hr - ft 



BTU 



2 



_ w 
t " AT 



783.39 



41 
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h r = <r * (T + ^ )(T 2 + T 2 ) 

= (0.1714 x 10" 8 ) (0.9) (571 + 530) (5?1 2 + 530 2 ) 



= 1 



BTU 



hr - ft 2 - °F 



h = 



T „ = 



h t - h r 



= 19.1 - 1 = 18.1 



BTU 



hr - ft - F 



T + T, 



<x> 



111 + 70 



= 90.5 F 



k„ = 0.0154 



BTU 



Nu = 



h *D 
c 



L f 



(evaluated at film temperature 

hr - ft - °F of 90.5 °F) 



= (18.1) (.. 4 I 

0.0154 



Nu = 294 at (9 = 63 degrees. 
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APPENDIX B 



UNCERTAINTY ANALYSIS 

Since there are uncertainties in the measured quantities, 
an uncertainty analysis was performed to ascertain the degree 
of uncertainty in the experimental results at three ranges of 
Reynolds numbers of 153.000, 90,000 and 16,000. The method 
of Kline and McClintock 6_7 was used. 

Being conservative, the following uncertainties were 



estimated in the measured variables: 

voltage, V £ 0.1 volts (20:1) 

resistance, R "t 0.4 ohms (20:1) 

temperature, T ± 1.0 °F (20:1) 

area, A t 0.001 ft. 2 (20:1) 

diameter, D t 0.001 ft (20:1) 

wall manometer, H t 0.05 in - H 2 o (20:1) 

radiation exchange 

factor - 0.05 (20:1) 



The uncertainties associated with the properties of air 
were considered to be negligible. 

The uncertainty analysis is performed for the data 
obtained at Reynolds numbers of 153.000. 90,000 and 16,000 
as follows; 

1 . Illustration of procedure of uncertainty analysis for 
Re = 153,000 : 
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Total Heat Transfer Coefficient, : 

As shown previously in Appendix A, the total heat 
transfer coefficient was calculated from the equation: 

, vf 

n t R • A • (T-Tco ) R . A • AT 



The uncertainty in the heat transfer coefficient is then 
expressed as: 






For the run at a Reynolds number of 153 » 000 , L/D = 8 , 
d/r = 4 and an angular location of 63 degrees: 

Y = 31 - 0.1 volts 

R =13.3-0.4 ohms 

AT = 41 t 1.0 °F ( T = 111°F = 571°R, ^ = 70 °F = 530°R) 
A = 0.314 t 0.001 ft 2 

H = 2.20 ± 0.05 in - H„o 



0 ), 



= 19.1 



BTU 



hr - ft 2 - °F 



(0 

h 

~h 



t 



0.04 or O) 

n t 



* 0,8 



BTU 

hr - ft 2 - °F 



This result does not take into consideration the 
uncertainty in the angular location of the liquid crystal 
band. This uncertainty is estimated to be - 3 degrees. 
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The uncertainty in the total heat transfer coefficient 
remained relatively constant around the circumference of the 
cylinder. However, the uncertainty in the angular location 
of the liquid crystal hand varied. An estimated maximum 
uncertainty of - 5 degrees occurred in the forward and rear 
stagnation regions and a minimum of - 1 degree occurred in the 
regions near the separation and reattachment points. 



Radiation Heat Transfer Coefficient, h r : 



The equation used to calculate h r was: 
h r = o’. (p [_ 2 -(T + T^ )(T 2 + t£ ) 



defining T m = ( T + T^ ) 



h r a * • "• • K -2 ■ T m 



Since a is constant, it is not considered in the uncertainty 
analysis. The uncertainty in h r is then expressed as: 




(jj. 



=0.06 or 



0)1 



= 0.06 



BTU 

hr - ft 2 - °F 



h = 1.00 t 0.O6 



BTU 



hr - ft* - °f 
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Local Nusselt number, Nuj 



The local Nusselt number v/as calculated from the 
equation : 



h -D 
Nu = 2 




Note that; h = h. - h 
c t r 



u) h 

~h 



c 



c 





Since the uncertainties in the properties of air are considered 
to be negligible and the uncertainty in the cylinder diameter 
is negligibly small; 




Nu = 294 - 13 
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Reynolds Number, Re : 



The equation used to calculate Reynolds number was: 



Re 



U -D 

CO 

u 



Therefore 

Re 
Re 

Since the uncertainty in the kinematic viscosity and in 
the cylinder diameter are negligible: 



Re U ro 

The quantity U^, is proportional to the square root of 
the wall manometer reading, hence: 

W Re _ 1 “k 




= 0.011 or m = 2000 

Re Re 

Re = 153.000 - 2000 

2 . The uncertainty analysis for the run at Re = 90,000: 
For the data obtained at a Reynolds number of 90 , 000 , 
L/D = 4, d/r = 2.0 and d = 67 degrees: 
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V 
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.0 


+ 


0 


.1 


volts 


13 


.0 


+ 


0 


.4 


ohms 


OO 


. 6 


+ 


1 


°F 


( T = 


0 . 


-t- 1 

-3- 

1 — i 


- 1 


0.001 


2 

ft 


0 . 


76 


+ 


0 


.05 


in - H 



-oo 



Similarly, the uncertainties calculated in following 
variables were: 



to. 



= 0.04 or co = ± 0.7 STIL 

h t h t hr - ft 2 - °F 



h t = 16.9 1 0.7 



BTU 



hr - ft 2 - °F 



to 



h 



, = 0 .06 or to, = 0 .06 

h r h r 



BTU 



hr - ft 2 - °F 



h r = 1.00 - 0.06 



BTU 



hr - ft 2 - °F 



to 



Nu 



Mu = °-° 4 ' + 0r “Nu = 11 



CO 



Re 

Re 



Nu = 259 T 11 
.032 



co Re = 3000 



Re = 90,000 t 3000 



532. 4°R) 
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3 • The un certainty analysis for the run at Re = 16,000 : 
For the data obtained at a Reynolds number of 16,000 
L/D = 4, d/r = 2.0 and © = 65 degrees: 

V =21.1 - 0.1 volts 

R =13.4- -0.4 ohms 

AT = 43 t 1 °F ( T= 111°F = 571°R , Tco = 68 °F = 5 28 °R) 
A = 0.314 - 0.001 ft 2 

The following uncertainties were calculated: 

U h 

h t = 0.039 or co = - 0.3 ITU 

t t hr - ft 2 - °F 

h. = 8.3 - 0.3 ITU 

t P n 

hr - ft^ - °F 



0)i 



-j—- = 0.06 or 10. = - 0.06 

r r 



BTU 



hr - ft 2 - °F 



h r = 1.00 t 0.06 



BTU 



hr - ft 2 - °F 



M Nu 

Nu 



0.04- or w Nu = - 5 
Nu = 120 - 5 



For this particular Reynolds number it was impossible to 
measure the free stream velocity using the wall manometer. 
The hot wire anemometer was therefore used to measure the 
velocity. Thus, the uncertainty estimated in was: 
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- 2 ft/sec 



( 20 : 1 ) 



(Ott 

^Re u °° 2 

Re U ro 10 



0.2 or w Re = 3000 



Re = 16,000 - 3000 



The high uncertainty in Uoo and consequently in the Reynolds 
number was primarily due to the unstable operation of the 
fan motor at its lowest operating point. 
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APPENDIX C 



WIND TUNNEL VELOCITY PROFILES 

For the three ranges of Reynolds numbers of 16,000, 

90,000 and 153,000 and with L/D = 4, velocity profiles were 
obtained in the tunnel test section (see Figure 18 through 
23). Data were collected at d/r = 0 and d/r = 5-33 (see next 
page for data) . 

The hot wire anemometer was used to probe the test 
section at various positive and negative Y - locations at 
90 degrees and 270 degrees, measured from the forward 
stagnation point on the cylinder. The velocity profile at 
the test section entrance was also measured at various 
locations and it was found to be fairly uniform. 

The turbulence intensity in the test section was obtained 
for a Reynolds number of 153 >000 and was found to be approxi- 
mately 0.2 fo. 
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DATA for the velocity profiles 
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SCALE: 1” = 3" 



Figure 18. Veloci ty profile at +h"p~+r>-r> +iT ~ \ ' v ^ ^ • * -■? 

J p top of the cylinder for d/r=0.0 
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SCALE: 1" = 3" 



•7 / / r- 7— -r /-■ J / r - y < z r 1 7 z — y z 7 7 ; 

Figure 19 . Velocity profiles at the top and bottom of the 
cylinder for d/r = 5*33* 
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SCALE: 1" = 3" 






Figure 20. ^locityprofil 
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SCALE: 1" = 3” 



Figure 21. Velocity profiles at the top and bottom of the 
cylinder for d/r = 5*33* 
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SCA1E : 1" = 3 " 



— /„ ' r ■ z: — ?-—r^r z. zr^z-^ — :/z z ry.-.- ~v — zr^zz zzz r ^: , — . — t rr? 
Figure 22. Velocity profile at the top of the cylinder for 
d/r = 0.0. 
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SCALE: 1" = 3" 



~ J—./ i • z • ■■■ z ir: -r t rr zLi'-jZii-.y- - .: ' 

Figure 23. Velocity profiles at the top and "bottom of the 

cylinder for d/r = 5«33< 
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APPENDIX D 



Additional plots of the Nusselt number versus gap size 
for Reynolds numbers of 16,000, 90,000 and 153 , 000 . 
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Figure 24. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r = 0.0. 
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Figure 25* Local Nusselt number on "the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r = 0.25- 
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Figure 2 6 . Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r =0.5. 
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Figure 27. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r = 1.0. 
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Figure 28. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r =2.0. 



79 



180 



Nu 




Figure 29. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r = 3.0. 
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Figure 30- Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r = 4.0. 
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Figure Jl. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 
Re = 16,000, L/D = 4, d/r = 5 -33. 
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Figure 32. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 90,000, L/D = 4, d/r = 0.0. 
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Figure 33* Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 90,000, L/D = 4, d/r = 0.25* 
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Figure 3^* Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 90,000, L/D = 4, d/r = 0.5« 
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Figure 35- Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 90 < 000 , L/D = 4, d/r = 1 . 0 . 
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Figure 36 . Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for -Re = 90,000, L/D = 4, d/r = 2.0. 
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Figure 37 * Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 90»000, L/D = 4 , d/r = 3*0* 
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Figure 38. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for He = 90 , 000 , L/D = 4, d/r = 4.0. 
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Figure 39- Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 90 , 000 , L/D = 4, d/r = 5*33* 
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Figure 40 . Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 > 000, L/D = 2, d/r = 0.0. 
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Figure 4l. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for He = 153 > 000 , L/D = 2, d/r = 0.25> 
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Figure ^2. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153.000. L/D = 2, d/r = 0.5- 
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Figure 43. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153.000, L/D = 2, d/r = 1.0. 
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Figure 44. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153.000, L/D = 2, d/r = 2.0. 
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Figure 4 5 . Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 > 000 , L/D = 2, d/r = 3*0. 
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Figure 46. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 .000, L/D = 2, d/r = 4.0. 
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Figure 47. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153>OOO f L/D = 2, d/r = 5*33* 
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Figure 48. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 > 000 , L /D - 8, d/r = 0.0. 
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Figure 49. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153.000, L/D =8, d/r = 0.25* 
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Figure 50. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 » 000, L/D = 8, d/r = 0.5* 
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Figure 51* Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153*000, L/D = 8, d/r = 1.0. 
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Figure 52. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 » 000 , L/D = 8, d/r = 2.0. 
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Figure 53- Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 >000, L/D = 8, d/r = 3*0. 
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Figure 54. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 > 000 , L/D = 8, d/r = 4.0. 
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Figure 55* Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 
surface for Re = 153 » 000 , L/D = 8, d/r = 5*33* 
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Figure 56 . Local Nusselt number on the surface of a free 
cylinder for Re = 16,000. 
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Figure 57 • Local Nusselt number and pressure coefficient on 
the surface of a free cylinder for Re = 90 , 000 . 
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Figure 58. Local Nusselt number and pressure coefficient 

on the surface of a free cylinder for Re = 153.000. 
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